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ABSTRACT 

Ultralight axions which couple sufficiently strongly to photons can leave imprints on the sky 
at diverse frequencies by mixing with cosmic light in the presence of background magnetic 
holds. We explore such direction dependent grey-body distortions of the CMB spectrum, 
enhanced by resonant conditions in the IGM plasma. We also find that if such axions are 
produced in the early universe and represent a subdominant dark radiation component today, 
they could convert into X-rays in supervoids, and brighten them at X-ray frequencies. 
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Axions are hypotetical pseudo-scalar particles motivated by - to date - the only known 
dynamical solution of the strong CP problem. In the Peccei-Quinn mechanism [I], an anoma¬ 
lous symmetry is spontaneously broken at a high scale. The phase of the complex held which 
picks a vev is the Goldstone mode, which gets a potential from nonperturbative QCD cor¬ 
rections, so that at the minimum the residual CP-violating phase is tiny, Ocp < 10- 9 PS], 
UV completions of the Standard Model, provided by string theory constructions, give rise 
to many more such phases, which at low energies behave as ultralight fields. One linear 
combination couples to QCD, and gets a mass from the QCD nonperturbative effects. Oth¬ 
ers may remain much lighter, and therefore active at cosmological scales. If they couple 
to long range 17(1) fields, in particular to electromagnetism, they could yield interesting 
astrophysical effects, providing extra dimming of supernovae, distorting CMB and affecting 
quasar spectra and cosmic rays - see for example m-ra- Even if they are decoupled from 
gauge theories, but only couple gravitationally, owing to their lightness they can still leave 
astrophysical signatures, in particular due to the phenomenon of superradiance as explained 
in [22][23]. 

Clearly, searching for axions is easier when there are direct axion-matter couplings. These 
couplings are given by irrelevant operators in order to maintain shift symmetry, and thus 
are suppressed by some high mass scale. However in the presence of nonvanishing back¬ 
ground fields these couplings can be enhanced by the background and yield possibly observ¬ 
able effects. For example, if one such ultralight axion couples to photons, by the standard 
dimension-5 operator L a _ 7 = , in the presence of background magnetic fields 

this will induce a bilinear term L a _ 7e g- = Bbacl ^ round (f>E that yields photon-axion oscillations. 
Depending on the coherence length of the background field, distance travelled and the envi¬ 
ronmental conditions, an incoming flux of photons can be partially converted to axions and 
evade detection. What remains, however, will - in general - have a different distribution 
over the range of frequencies than it had at the source, since the mixing rate is frequency 
dependent. Hence if a source of photons were well understood, seeking for the spectral dis¬ 
tortions can either indirectly reveal the presence of a photon-axion mixing, or constrain it. 
Conversely, it is also possible that in the presence of relic axions, the conversion can in fact 
enhance a flux of photons, making certain sections of the sky where the mixing is strong 
brighter than the others. From a fixed observer’s perspective, this would make the spectral 
distortions line-of-sight dependent, and so anisotropic. 

These phenomena were applied to the search for ultralight axions, with masses < lCC 15 eV. 
The reason is that in cosmological environments, the presence of the plasma, with plasma 
frequency induced effective photon mass, m 7 > 10~ 15 eV, very efficiently screens the mixing. 
Similarly a large axion mass would screen it even more efficiently. However these screening 
effects may cancel, enhancing the mixing to close to maximal even for low frequencies. This 
could happen if the plasma frequency is sufficiently close to the axion mass in large regions 
of space, where the plasma is sufficiently smooth. Candidates for such regions, which should 
be quiet, are very large voids. In light of the fact that such supervoids have been recently 
discovered [23} [25], we will consider what happens with the photons and axions inside them, 
if the magnetic fields are large and the plasma sufficiently slowly varying. Also, we will 
consider the mixing effects in the clusters, where the plasma frequency is higher than in the 
voids, but the magnetic fields can also be much larger. As long as the magnetic coherence 
length is comparable to the length scale controlling the variation of the plasma density, the 
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conditions for resonant mixing may be met. For this reason we will consider axions with 
masses in the range 10~ 12 — ICR 16 eV, with the mass being a scanning parameter, such that 
the resonant conditions may be met somewhere in space. This can occur in the framework 
of the axiverse, with a plenitude of light axions [22J. Of course, it is still necessary that the 
axion which can be in resonance is the one coupling to the usual electromagnetism. 

Assuming this occurs, we will show how the mixing can produce a direction-dependent 
grey-body distortion of the CMB spectrum. While our results are to be viewed as estimates, 
they are an interesting indicator of the potential of cosmological observations for discovering 
axions in this as yet largely unexplored regime of parameters, where the observations are 
rapidly improving. So, for example, a dedicated analysis on the Planck data can provide a 
map of the regions through which photons have disappeared into axion^J 

Second, we will consider the possibility that there is a residual high frequency axion 
background, as suggested in [26]. The premise of this scenario is that the universe contains 
axion ‘dark radiation’, generated by a late decay of heavy string moduli, with a spectrum 
peaked today at the soft X-ray frequencies. These axions would convert into X-ray photons 
in background magnetic fields inside clusters. This yields an excess X-ray luminosity from 
galaxy clusters 123 1231 E9], claimed to match some recent observations [30]. While the 
uncertainties in the predicted X-ray luminosity at present limit the direct discovery potential 
of such searches, we point out that a possible tie-breaker could be provided by voids. Namely, 
inside the magnetic fields in voids, which are weaker but possibly more coherent, the relic 
axions could convert relatively efficiently into photons, making the void brighter in the 
X-ray part of the spectrum. Such a signal might be present in the Rosat X-ray and XMM- 
NEWTON data, and a dedicated analysis therefore seems warranted. 

The axion-photon mixing dynamics is described by the action 


S 





( 1 ) 


where is the standard electromagnetic tensor, F^ u = 2 J-g elWpa Fpm m a is the axion mass 
and M is a mass parameter controlling the coupling strength. We will be picking M near 
the current limits, M ~ 10 11 GeV [BT], since our main goal is to illustrate the sensitivity 
of the current observations to the effects induced by ultralight axions. So we assume that 
the couplings and environmental magnetic fields are as strong as possible while not being 
excluded by observations at present. The future data may constrain these couplings more 
tightly. We ignore the space-time curvature effects because we will work with the axion 
masses m a much greater than the relevant curvature scales. Furthermore, the flavor mixing 
can only take place in relatively low plasma density environments with significant background 
magnetic fields, at least of the order of nano-Gauss, but with coherent lengths much shorter 
than the curvature scales. Such conditions might only occur at very low redshifts < 2, and 
so the cosmic evolution is largely negligible. This also implies, that any effects from photon- 
axion mixing will be suppressed for very old light, since at such early times the plasma 
frequencies were higher and magnetic fields likely negligibly small according to current ideas 
on magnetogenesis. 

1 We are aware that this analysis may be difficult because of calibration uncertainties. However, future 
experiments dedicated to measuring the CMB spectrum will be sensitive to the photon disappearance. 
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In the presence of a background magnetic field B, the equations yield a simple bilinear 
system which describes photon-axion oscillations 013 E], The conversion probability of a 
beam of photons which travels a distance z in the presence of a magnetic held B admits a 
perturbative solution 

2 

( 2 ) 

where A a = m 2 /(2u), A p = uj 2 / 2u, (i\\ = B\\/M, and oj 2 = A7ian e /m e is the plasma 
frequency, with n e the number density of free electrons. In the integral above we allow for 
spatial variations of the magnetic held B and the electron number density n e , but treat it 
adiabatically as in [6], 

To gain physical insight into the dynamics, we can specialize to the case of constant 
coherent magnetic held filling a spatial domain of size L do m and containing a constant electron 
number density n e . Solving the coupled system of held equations governing photon-axion 
propagation [7j, we find the standard transition probability of a 2-level quantum system: 

D _ h 2 Llom ( sin(r ]L dom )\ 2 ^ 

* - ~r“ l. nL^ ) ’ (3) 



./ — 7772 ^2_|_^^2^j2 

where rj is the inverse of the oscillation length, rj = ——-—^ -. Depending on the 

parameters and the frequency of the modes, several different regimes can set in: 


\(jj2 _777,^ I 

1. maximal mixing , with oj p 0fl ° everywhere in the space where photons propagate, 

and the probability of transitions P t — sin 2 (/iL ( iom/2) is independent of frequency 
regardless of the size of magnetic domains; 


2. non-maximal mixing, where gL^om 1, which means that for the hxed frequency of 
interest the magnetic domain is very thin: in this case, Pt — ia 2 L 2 dom /A is constant and 
frequency- independent; 


3. non-maximal mixing, with T]L dom 1, implying that for a hxed frequency the magnetic 
domain is very thick: in this case, the oscillations can be averaged and we end up with 
Pt — /j 2 /(8r] 2 ) ~ 2h 2 oj 2 /(uj 2 — m 2 ) 2 , which has a co 2 frequency dependence, and a 
small amplitude; the crossover frequency between the latter two regimes is 0J C rossover — 
— ™ a)L * , where L* = min (L dom ,M/B) [21]; 

4. resonant regime, where in some regions | oj 2 — m 2 \ happens to be very small compared 
to 2 (loj locally, due to the fact that the plasma density changes to be close to the axion 
mass. In this case, inside this region the conditions for maximal mixing are fulfilled and 
frequency dependence is very weak locally; the photon survival probability inside only 
such regions is P t ~ sin 2 (/iL dom /2), while outside is given by the formulas governing 
case 2) or 3) 0 152] . 


So axions will mix with photons of any frequency at some level in background magnetic 
fields, with the transition rate depending on the frequency. This has possible important 
cosmological implications for the Cosmic Microwave Background photons (CMB) and for 
the high energy cosmic ray photons. We will explore them below, focusing in fact on the 
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case 4) for both CMB photons and soft X-rays in clusters, and case 1) for high energy X- 
rays. Note, of course, that in a realistic situation the actually observed flux of photons will 
be the result of propagation through many different domains along the line of sight. This 
means that the total rate of photon-axion conversions will be determined by a combination 
of different dynamical regimes, possibly involving - in general - a combination of some, or 
even all, of the cases listed above, due to the random variation of environmental parameters, 
including strength and direction of background magnetic fields, plasma density, residual 
axion abundance and so on. For example, the photon-axion conversion in the optical range, 
which can affect super novae, is a consequence of integrating over many randomly oriented 
domains where some of the conditions listed above apply, as described in jTj. For CMB, and 
in the environments of interest to us in what follows, however, the transitions in individual 
domains are suppressed unless the special conditions allowing faster transition rates occur. 
Hence we will ignore the contributions from many domain in the rest of this article, focusing 
exclusively on the cases where one domain along the line of path contributes dominantly to 
photon-axion conversion. 

CMB has a thermal spectrum to a very good accuracy because of tight coupling with 
baryons before recombination. This form of the distribution is preserved by the adiabatic 
expansion of the Universe. The temperature varies very little with direction, because inho¬ 
mogeneities were small at recombination. Perturbations in the local density of photons can 
be interpreted as a local variation in the black-body temperature: 


Ibb (w, h) 


2 hu 3 
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where Ibb( w,h) is the radiation spectrum as a function of energy E — huo and angular 
direction n. If there is any local variation of the temperature, there will be a local excess or 
lack of photon energy density, given by the leading order variation of Q 


51 (u>, h) 


dl 5T(n) x 5T(n) 

dWT Tq ~J\T ~ BB l-e~* T 0 
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where x = huj/(kBT 0 ) and T 0 = 2.726 K is the direction-averaged temperature measured by 
COBE-FIRAS. Since in the presence of background magnetic fields and axions CMB photons 
can disappear, at a rate which depends in general on their frequency, and varies with the 
line of sight, local spectral distortions can be induced. Indeed, the spectrum of photons 
surviving the propagation from the last scattering surface to the detector is 


I c (u>, ft) = I BB (u)(l - P c (n,u)), 


( 6 ) 


which obviously deviates from a pure black body spectrum. Converting this spectral distor¬ 
tion into thermodynamic temperature units, we get 

5T 1 - e~ x 

= —P c (n,Lo) -. 7 

1 x 

This shows that the thermodynamic temperature acquires a frequency dependence even if 
the conversion probability is independent of frequency - as in the case 4) listed above, which 
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Figure 1: Plot of the function (1 — e x )/x , as a function of the frequency in GHz. The points 
correspond to the central values of the Planck frequency bands. 


we will focus on since if it is realized for the CMB photons and realistic weak intergalactic 
magnetic fields, it yields the dominant distortions of the spectral distribution. The cases 1) 
and 3) can never yield significant effects for CMB photons with realistic parameters taken 
into account. The case 2) is more interesting, but also generically small without resonance 
conditions (however see [21j for possible thermal distortions in extragalactic space which are 
potentially discernible by comparison of several different frequency channel data sets).The 
resultant frequency dependence of the distorted spectrum due to this effect is shown in Fig. [lj 

Let us now demonstrate that the scale of the effect is within the reach of observations, 
as claimed above. The weakest assumption in a part of our argument is that the magnetic 
fields should be sufficiently strong. Since they are not known directly we will assume that 
they are close to saturating the current limits. A justification of this assumption is that 
while the size of the effects we are discussing might not be viewed as a firm prediction, it 
is a qualitative prediction of a type of an imprint in the CMB. If observed, such effects 
could be independent checks not only of the axion’s presence but also of the extragalactic 
magnetic fields, which are otherwise very hard to observe directly. Alternatively, the effect 
could also be enhanced even if the magnetic fields are weaker, if the resonance conditions 
(case 4) above are met in some sufficiently large and smooth region of space with coherent 
magnetic fields. This at first sounds like fine tuning of the Lagrangian parameters - i.e., the 
axion mass m a - to the environmental parameters - i.e., the plasma frequency oj p (10). One 
may worry if this can ever be achieved. A viable way to relax this criticism is the paradigm 
of the axiverse [22] which posits that there exist many light axions, with masses that are 
small and that scan the full range of cosmologically interesting scales, down to the current 
value of the Hubble parameter Ho- In such a framework, the resonant conditions are realized 
along some directions for at least some of the axions. The signatures which we discuss can 
arise if that particular axion is the one which couples to the usual electromagnetism. 

With this in mind, we first consider the photon-axion conversion in clusters. There occur 
the coherent magnetic fields B ~ /rG, with coherence length ~ kpc. This implies that the // 
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The plasma frequency can be estimated as 


UJp 


= 5.5 x 1CT‘V 2 (A) eV , 


( 9 ) 


where x e is the ionization fraction of electrons (x e ~ 1 at low redshifts), and A = p m /p is the 
local overdensity parameter. For a cluster, we expect A ~ 10 3 — 10 4 . The angular frequency 
for a typical microwave photon of frequency z/ 100 = z//100GHz is oj = 4.1 x 10~ 4 z/ 100 eV. 
Using these values, and defining B x = B\\/X, M n = M/10 41 GeV, the (square of the) 
inverse oscillation length will be (with axion mass nn a —> m a /eV in what follows) 
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Since the product r]L dom is generically large, the photon-axion mixing is very suppressed 
away from the resonance. At the resonance, however, the transition probability is 

P t , res = sin 2 (^L dom ) -2.9x10 3 , ( n ) 

which in fact is even larger than the observed CMB anisotropies! Clearly, this sounds like 
a dangerously large effect. But with slightly weaker and less coherent magnetic fields, with 
weaker couplings, and less uniform plasma the conversion rates are suppressed and the 
signatures are weaker. In fact one would expect this anyway, because the variation of the 
plasma density spoils the resonant conditions. Also, the conversion factor to thermodynamic 
temperature is always less than unity (Fig. [Tj). 

We note that there is another interesting possibility where the photon-axion transitions 
can occur through voids, which are actually the largest fraction of the volume of the Universe. 
This can happen if the magnetic fields are large enough; assuming the values B ~ nG, 
L dom ~ Mpc and an underdensity A ~ 10 _1 , we get 
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As in the case of clusters, r)L dom 1, and again the effect should happen in resonance in 
order to maximize it. Since the plasma density is about co p i = 1.2 x 10~ 14 y /f pp||^=(jeV, the 
resonant regime can occur for light axions, with masses m a A 10~ 14 eV. While the required 
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resonance regime is narrower in the parameter space, it can nevertheless occur in random 
voids. The size of the transition probability can actually be of the same order as for clusters: 

= 2.9xl0-(^) 2 (^) 2 . (13) 

So as a consequence, the thermal distortions could be imprinted as a signature of voids 
bearing large magnetic fields. Again, the scale of the transition probability is larger than the 
observed CMB anisotropies, but that is a consequence of our assumptions about the strength 
of photon-axion couplings, magnetic field strengths and coherence, and plasma homogeneity 
scales. In reality, these can be all weakened yielding smaller but possibly observable effects. 

Finally, we stress that to be observable for realistic parameters, the spectral distortions 
from the photon-axion conversions should be distinguishable from the distortions due to 
the foregrounds. The foreground effects are typically expressed in the units of antenna (or 
brightness) temperature, defined in the Rayleigh-Jeans approximation, which yields 

TA(u ' v) = 2t^ I(u ' v) - (14) 

In the antenna temperature units, the distortion due to axion conversion translates into 

— = -P c (n,tn) —(15) 

A measurement of such spectral distortions may be feasible with the actual Planck data (33] . 
and would provide a map of photon “disappearance” along the line of sight. The possible 
degeneracy with foregrounds, in particular synchrotron at low frequencies, may be resolved 
by noting that the spectral index of the axion-induced corrections f3 ~ —0.3 around 30 GHz 
is quite different from the synchrotron distortions, and so the foreground effects may be 
subtracted away at least in principle. 

Similar phenomena also arise for very high energy photons. These photons - i.e. cosmic 
X-rays and y-rays - could be less sensitive to electron plasma since their frequency is so 
high. In extragalactic space, they are in the regime of maximal mixing, i. e. their dynamics 
falls in the category of the case 1) we listed above. In clusters, where the plasma density 
is higher and magnetic fields are stronger, the transition between the maximal and non 
maximal mixing can occur at higher frequency. This is further exacerbated if the photon- 
axion coupling is weaker as considered by [26]. However, we note that if the axion mass is 
in the resonant regime relative to the plasma, as described by the case 4) which we listed 
above, the transitions can be enhanced yielding a stronger signal. As we noted above, the 
possibility of photon-axion resonance in plasma can be naturally realized in the context of 
the axiverse where the axion mass parameter scans a range of scales, and so it could hit the 
right value to be in resonance in some regions of space. This makes the X-ray signatures in 
both extragalactic space and in clusters possibly sensitive to the presence of such axions. 

In more detail, the enhanced X-ray luminosity can occur if there is a relic cosmic axion 
background (CAB), of very high energy axions. Specifically, [26] proposed a scenario where 
such relic background axions arise as a subdominant dark radiation from a late decay of 
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heavy string moduli dark matter, 
order 

Pcab 


The estimated density of CAB axions in [26J is of the 


= 1.6 x 10 60 


erg A N eff 
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where AN e jf is the fractional contribution of the axions to the effective number of relativistic 
species in addition to photons and Standard Model neutrinos. Such radiation amounts to 
only 10% of the CMB energy density, however with the peak around keV energies. As argued 
in [27], in the presence of a background magnetic field these axions will convert into soft X- 
ray photons. The expression for conversion probability is of course eq. ([3]). The parameters 
which the authors of ra employ in their analysis correspond to the dynamics of case 2), 
with marginally non-maximal mixing and frequency sensitivity. 

However, as we noted above, if the resonant regime occurs for the conversion of CMB 
photons, it can also occur for the conversion into X-ray photons. Obviously, this will not be 
realized inside every cluster, or over the whole volume of a cluster where the resonance may 
take place. Nevertheless, since the environmental conditions vary from cluster to cluster, 
and the axion mass can scan the regime of scales, the resonant conditions can be satisfied 
somewhere. The conversion rate per axion per second will be r ~ P C onversion c/ Ld om - Using 
our previous estimates of the photon-axion conversion probabilities, we can estimate the 
total X-ray luminosity density of a cluster where resonant conversion sets, 
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The overall scale of the excess luminosity in (17) by itself is too large . However, as be¬ 


fore with the imprints on CMB, it diminishes with weaker photon-axion couplings and less 
coherent plasma and magnetic fields. Specifically, since the total X-ray luminosity of a typ¬ 
ical cluster due to conventional electromagnetic processes is L ~ 10 44 ergs -1 , taking the 
weaker photon-axion couplings Mu ~ 0(1000), and without assuming resonance, the au¬ 
thors of [27] conclude that conversion of the CAB axions into photons can yield an observable 
X-ray luminosity excess if the relic axion abundance is sufficiently high, and even point to 
certain reported astrophysical anomalies. However, because of astrophysical uncertainties, 
it is harder to firmly establish that the axions are the cause for the excess luminosity noted 
in those cases [2Ej. Nevertheless, the possibility remains, and as we see above, the resonance 
conditions - if realized - can yield effects with even weaker couplings. 

Here we also point out another intriguing possibility for the detection of high energy 
CAB axions: instead of looking at the clusters, look at the voids. The point is that although 
in the voids the magnetic field is weaker, the axion density is parametrically the same, since 
they are ultra-relativistic. With nano-Gauss magnetic fields in large voids with the resonant 
conditions satisfied in regions of size Ld om ~ Mpc, we can estimate the luminosity by scaling 
down the formula ( [l7| , using the fact that the magnetic fields are a thousand times weaker 
while the domain is a thousand times larger. This yields 
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which is also very large. Again, as we stressed above, this is not a problem of principle, 
because the luminosity drops when the magnetic fields are weaker, the resonant domains are 













smaller than the whole ~ Mpc sized underdensity, and the axion relic abundance is smaller. 
Nevertheless our point that this is a possible signature of the photon-axion mixing in the 
voids remains. This could be a new portal into the dark sector of the universe, and deserves 
more attention, both theoretical and observational. 

Yet another interesting signature of photon-axion conversion, and evidence of axions, 
is at yet higher energies, for 7 -rays. These can be generated in high energy processes in 
astrophysical sources such as AGNs. In the standard picture, the spectrum of 7 -rays drops 
sharply at energies > ITeV because of pair production in collisions with extra-galactic 
background light. Previously, the photon-axion mixing mechanism at such high frequency 
has been considered to allow ultra-high cosmic 7 -rays from very distant astrophysical sources 
to surpass the GZK cutoff limits ra- However, as suggested in pans], photon-axion mixing 
could harden the observed high-energy spectrum. Indeed, as explained in [13], high energy 7 - 
rays will partly convert into axions in the strong magnetic held close to the source, and travel 
undisturbed to our Galaxy, where they will convert back to photons and reach our detectors. 
While the astrophysical uncertainties obscure the precise identification of such events, it may 
be possible to get to them statistically, by systematically searching for anisotropies in the 
spectra of AGN’s. In particular, [T9] proposes exploring correlation between sources on small 
scales and an anti-correlation on large scales, which is a common characteristic of different 
magnetic held models. While the tests carried out on current Fermi data are inconclusive at 
the moment, the future observations by the Cherenkov Telescope Array may have enough 
statistics to test part of the parameter space of photon-axion mixing and the presence of 
light axions which couple to photons. 

To summarize, we have considered observable effects of photon-axion mixing in cosmo¬ 
logical conditions. We have argued that the photon-axion mixing can leave an imprint in 
the CMB, especially for the case where the axion mass - as in the axiverse - could happen 
to have the right value in order to be in resonance with plasma frequency in voids. Future 
explorations of CMB may be sensitive to such imprints, as already argued in ED. where the 
case was made that transitions out of resonance could yield distortions of the CMB. Here 
we pointed out that distortions of the overall temperature arise even when the dominant 
photon-axion conversions are frequency independent, just because the survival probability 
of a CMB photon is changing the normalization of the total black body radiance. Such 
effects may be the strongest in the voids, and given that now very large voids appear to be 
an experimental fact [22 [25], it seems warranted to look for such effects. It may be possible 
to determine such temperature shifts by the frequency-dependent analysis of the observed 
spectra. We also noted that the voids could be glowing in the X-ray spectrum if there 
is a relic cosmic axion background as well. Given the difficulties in directly searching for 
weakly coupled dark sectors such as axions, these cosmological windows may be a promising 
direction to explore seriously. 
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